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Abstract 

In the recently proposed minimal supersymmetric SO(10) model, the neutrino Dirac 
Yukawa coupling matrix, together with all the other fermion mass matrices, is completely 
determined once free parameters in the model are appropriately fixed so as to accommo- 
date the recent neutrino oscillation data. Using this unambiguous neutrino Dirac Yukawa 
couplings, we calculate the lepton flavor violating (LFV) processes and the muon g — 2 
assuming the minimal supergravity scenario. The resultant rates of the LFV processes are 
found to be large enough to well exceed the proposed future experimental bound, while 
the magnitude of the muon g — 2 can be within the recent result by Brookhaven E821 
experiment. Furthermore, we find that there exists a parameter region which can simul- 
taneously realize the neutralino cold dark matter abundance consistent with the recent 
WMAP data. 
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1 Introduction 



The problems of the observed solar neutrino deficit and the atmospheric neutrino anomaly 
can be naturally solved through the interpretation of the neutrino oscillations among the 
different flavor neutrinos. Now the evidence of the neutrino masses and flavor mixings are 
almost established [I], and this is also the evidence of new physics beyond the standard 
model. Interestingly, the neutrino mass and mixing properties have been revealed to be 
very different from those of the other fermions, namely, neutrino masses are very small 
and the flavor mixing angles are very large. A new physics must explain them. 

The supersymmetric (SUSY) SO(10) grand unified model is one of the attractive 
candidates for new physics. The experimental data support the unification of the three 
gauge couplings at the grand unified theory (GUT) scale M G ~ 2 x 10 16 GeV with the 
particle contents of the minimal supersymmetric standard model (MSSM) [2 J. This model 
incorporates the right-handed neutrinos as the member of the 16 representation together 
with other standard model fermions, and provides the natural explanation of the smallness 
of the neutrino masses through the see-saw mechanism 

However, it is a non-trivial problem whether a SO(10) model can simultaneously ac- 
commodate all the observed quark-lepton mass matrix data, since the quark and lepton 
mass matrices are related with each other and severely constrained by virtue of the quark- 
lepton unification into the 16 representation. There are lots of works on this subject. 
Recently it has been shown [3] that the so-called minimal supersymmetric SO (10) model 
can simultaneously reproduce all the observed quark-lepton mass matrix data involving 
the neutrino oscillation data by fixing free parameters in the model appropriately. In this 
model, the neutrino Dirac Yukawa coupling matrix, together with all the other fermion 
mass matrices, is completely determined. 

The evidence of the neutrino flavor mixing implies that the lepton flavor of each 
generation is not individually conserved. Therefore the lepton flavor violating (LFV) 
processes in the charged-lepton sector such as fi — > ej, r — > /ry are allowed. In simply 
extended models so as to incorporate massive neutrinos into the standard model, the rate 
of the LFV processes is accompanied by a highly suppression factor, the ratio of neutrino 
mass to the weak boson mass, because of the GIM mechanism, and is far out of the reach 
of the experimental detection. However, in supersymmetric models, the situation is quite 
different. In this case, soft SUSY breaking parameters can be new LFV sources, and the 
rate of the LFV processes are suppressed by only the scale of the soft SUSY breaking 
parameters which is assumed to be the electroweak scale. Thus the huge enhancement 
occurs compared to the previous case. In fact, the LFV processes can be one of the most 
important processes as the low-energy SUSY search [3]. 

Any mechanism of the SUSY breaking mediation are normally considered so as to 
induce the soft SUSY breaking parameters being almost flavor blind and real because of 
sever constraints on the SUSY flavor and CP violating processes jHj- Therefore the rate 
of the LFV processes is negligible at the SUSY breaking mediation scale. However, in a 
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model with the LFV interactions, the sizable LFV sources at low energies may be induced 
through the renormalization effects. As well-motivated candidates for such models, the 
SUSY GUTs [Zj and the models with see-saw mechanism jH] have been considered. 

The magnitude of the LFV rate depends on the SUSY breaking mediation scenario 
and the LFV interactions. In this paper, we evaluate the rate of the LFV processes in 
the minimal SUSY SO(10) model, where the neutrino Dirac Yukawa couplings are the 
primary LFV sources. The minimal supergravity (mSUGRA) scenario is assumed as the 
SUSY breaking mediation mechanism. As mentioned above, all the Yukawa couplings are 
determined in the minimal SUSY SO(10) model, we can give the definite predictions for 
the rate of the LFV processes. 

Recent result of the muon g — 2 from Brookhaven E821 experiment |3J reported the 
deviation from the standard model prediction. Although the uncertainty of the standard 
model prediction has not been settled ^U], we would expect that the deviation lies in the 
range 

Sa^ = a M (£821) - a^(SM) = (10 - 40) x 10" 10 , (1) 

where a M = (g M — 2)/2. This may be the first evidence of the low-energy SUSY [11. 
Except for the flavor indeces, both of the LFV processes such as fi — > and the muon 
g — 2 originate from the same dipole- moment operators, and there is a correlation between 
their magnitudes [T2]. In the following, we calculate the SUSY contribution to the muon 
g — 2 with the same input parameters as in the calculations of the LFV processes, and 
show their correlations. 

This paper is organized as follows: in the next section, we give a brief review of the 
minimal SUSY SO(10) model developed in [3], and present the predicted neutrino Dirac 
Yukawa coupling matrix which plays the central role in our results. In section 3, using 
the neutrino Yukawa coupling matrix, we calculate the rate of LFV processes and the 
magnitude of the SUSY contribution to the muon g — 2 for various input parameters in 
mSUGRA scenario. The last section is devoted to summary. 

2 Minimal SO (10) model and its predictions 

Here we give a brief review of the minimal SUSY SO(10) model 1 recently reconsidered in 
Ref. [3], and summarize its predictions. Even when we concentrate our discussion on the 
issue how to reproduce the realistic fermion mass matrices in the SO(10) model, there 
are lots of possibilities for introduction of Higgs multiplets. The minimal supersymmetric 

1 There are other models which are also called "minimal SUSY SO(10) model" in literatures 
Such models include only 10 Higgs multiplet and realize the Yukawa unification between bottom quark 
and tau lepton. In our "minimal" model, a 126 Higgs multiplet is additionally introduced, and Yukawa 
unification is not necessarily satisfied. 
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SO(10) model is the one where only one 10 and one 126 Higgs multiplets have Yukawa 
couplings (superpotential) with 16 matter multiplets such as 

W Y = Yiil6iH 10 16j + F^ 6 16^ 126 16,- , (2) 

where 16, is the matter multiplet of the i-th generation, H w and Hi 2 § are the Higgs 
multiplet of 10 and 126 representations under SO(10), respectively. Note that, by virtue 
of the gauge symmetry, the Yukawa couplings, Y\q and Yi 2 q, are, in general, complex 
symmetric 3x3 matrices. 

We assume some appropriate Higgs multiplets, whose vacuum expectation values 
(VEVs) correctly break the SO(10) GUT gauge symmetry into the standard model one 
at the GUT scale, Mq ~ 2 x 10 16 GeV. Suppose the Pati-Salam subgroup [Hj, G422 = 
<S77(4) C x SU(2)l x SU(2)r, at the intermediate breaking stage. Under this symme- 
try, the above Higgs multiplets are decomposed as 10 — > (6, 1, 1) + (1, 2, 2) and 126 — > 
(6, 1,1) + (10, 3,1) + (TO, 1,3) + (15, 2, 2), while 16 -> (4,2, 1) + (2, 1,2). Breaking down 
to the standard model gauge group, SU(4) C x SU (2) R — > SU (3) c x U(1)y, is accomplished 
by non-zero VEV of the (10, 1,3) Higgs multiplet. Note that the Majorana masses for 
the right-handed neutrinos are also generated by this VEV through the Yukawa coupling 
F 126 in Eq. 0. 

After this symmetry breaking, we find two pair of Higgs doublets in the same repre- 
sentation as the pair in the MSSM. One pair comes from (1,2,2) C 10 and the other 
comes from (15, 2, 2) C 126. Using these two pairs of the Higgs doublets, the Yukawa 
couplings of Eq. (J2J) are rewritten as 

W Y = (t£), (YftHZ + Qj + (c$), + Y&H*n) q 5 

+ ( u R)i (^10-^10 ~~ SYjgif " 26 ) tj + (e c R )i (Y$Hi — SY^H^^j £j 
+ (^)i (Y& v R ) , (3) 

where ur, cIr, Vr and e R are the right-handed SU(2)l singlet quark and lepton super- 
fields, q and £ are the left-handed SU(2)l doublet quark and lepton superfields, Hiq and 
H126 are up-type and down- type Higgs doublet superfields originated from H w and H 12 q, 
respectively, and the last term is the Majorana mass term of the right-handed neutrinos 
developed by the VEV of the (10, 1, 3) Higgs, v r. The factor —3 in the lepton sector is 
the Clebsch-Gordan coefficient. 

In order to preserve the successful gauge coupling unification, suppose that one pair of 
Higgs doublets given by a linear combination and H^i is light while the other pair 
is heavy (> Mq)- The light Higgs doublets are identified as the MSSM Higgs doublets 
(H u and Hd) and given by 

H u = a u H^Q + /3 u H^2q 

H d = a d H d w +p d H d 12&1 (4) 
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where a Ujd and (3 U}d denote elements of the unitary matrix which rotate the flavor basis 
in the original model into the (SUSY) mass eigenstates. Omitting the heavy Higgs mass 
eigenstates, the low energy superpotential is described by only the light Higgs doublets 
H u and H d such that 

W Y = (t£)i (a u Y^ + H u q 3 + (d%) t (a d Y^ + (3 d Y^ 6 ) H d q 3 

+ (yfti (a u YH - 3/9"y$,) HJ 3 + (e^ (a d Y^ - 3(3 d Y^) H d £ 3 

+ {y&Pr) , (5) 

where the formulas of the inverse unitary transformation of Eq. (J3J, = a u ' d H Utd + ■ ■ ■ 
and H12Q = (3 u,d H U)d + ■ • •, have been used. Note that the elements of the unitary matrix, 
a u,d and /3 u,d , are in general complex parameters, through which CP-violating phases are 
introduced into the fermion mass matrices. 

Providing the Higgs VEVs, H u = vsm/3 and H d = vcos/3 with v = 174GeV, the 
quark and lepton mass matrices can be read off as 2 



M u 


= C10M10 + c 12 eM 12e 


M d 


= M 10 + M 126 


M D 


= C W M W - 3Ci2 6 Mi2 6 


M e 


= M 10 -3M 126 


M R 


= C R M 12 6 , 



where M u , M d , M D} M e , and M R denote the up-type quark, down-type quark, neutrino 
Dirac, charged-lepton, and right-handed neutrino Majorana mass matrices, respectively. 
Note that all the quark and lepton mass matrices are characterized by only two basic mass 
matrices, M w and Mi 2 6, and three complex coefficients cio, cvi% and cr, which are defined 
as Mio = Y w a d v cos p, M 126 = Y 126 f3 d v cos (3, c w = (a u / a d ) t&n (3 , c 126 = (p u /p d ) tan (3 
and cr = VR/(j3 d v cos/3)), respectively. Except for cr, which is used to determine the 
overall neutrino mass scale, this system has fourteen free parameters in total ^] , and the 
strong predictability to the fermion mass matrices. 

Thirteen electroweak data of six quark masses, three mixing angles and one phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and three charged-lepton masses are 
extrapolated to the GUT scale according to the renormalization group equations (RGEs) 
with given tan/3, and are used as inputs at the GUT scale. 3 Solving the GUT mass 

2 In general, the SU{2)l triplet Higgs in (10,3,1) C 126 can obtain the VEV induced through the 
electroweak symmetry breaking and may play a crucial role of the light Majorana neutrino mass matrix. 
This possibility, so-called type II see-saw model, has been discussed in detail in Ref. ^3 ■ 

3 In our analysis, we have neglected SUSY threshold corrections to the quark Yukawa coupling whose 
importance in the case of large tan/3 has been pointed out ^7j. In our following analysis with tan/? = 45, 
we find the corrections about 10%, which would affect our final results. However we can check that this 
effects can be compensated away by slight changes of other input values, for example, mixing angles in 
CKM matrix within the experimental errors, and our final results is found to be almost unchanged. 
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matrix relation among quarks and charged-leptons obtained by Eq. (jUJ), we can describe 
the neutrino Dirac mass matrix Mb and M 12 6 as functions of only one free parameter a, 
the phase of the combination (3io + ci26)/(— cio + ci26)- The light Majorana neutrino mass 
matrix, M v = -M^M^M D = -c^M^M^Mn, given through the see-saw mechanism 
is extrapolated down to the electroweak scale according to the RGE for the effective 
dimension-five operator and is compared to the neutrino oscillation data for various 
cr. It has been shown [3] that the appropriate value of o can reproduce three neutrino 
flavor mixing angles in the Maki-Nakagawa-Sakata (MNS) matrix and the ratio of mass- 
squared differences ( Am|,/Am^) consistent with the current neutrino oscillation data. 
Here Am^ and Am^ are the mass-squared differences relevant for the solar and the 
atmospheric neutrino oscillations, respectively. The parameter cr is used to determine 
the overall scale of the right-handed Majorana neutrino mass, and thus the scale of Am^. 
Once the appropriate values of a and cr are chosen, all the fermion mass matrix (or 
fermion Yukawa couplings) are completely determined. Therefore the model has definite 
predictions to physics related to the fermion Yukawa couplings. 

Although in Ref. [3] various cases with given tan f3 = 40 — 55 have been analyzed, we 
consider only the case tan/? = 45 in the following. Our final result in the next section 
is almost insensitive to tan/? values in the above range. The predictions of the minimal 
SUSY SO (10) model that we need in our calculation in the next section are as follows 
j3]: with a = 3.198 fixed, the right-handed Majorana neutrino mass eigenvalues are found 
to be (in GeV) M Rl = 1.64 x 10 11 , Mr 2 = 2.50 x 10 12 and M Ri = 8.22 x 10 12 , where 
cr is fixed so that Am^ = 2 x 10~ 3 eV 2 . In the basis where both of the charged-lepton 
and right-handed Majorana neutrino mass matrices are diagonal with real and positive 
eigenvalues, the neutrino Dirac Yukawa coupling matrix at the GUT scale is found to be 



/ -0.000135 - 0.00273i 0.00113 + 0.0136i 0.0339 + 0.0580* 

0.00759 + 0.0119* -0.0270 - 0.00419i -0.272 - 0.175i j . (7) 
V -0.0280 + 0.00397i 0.0635 - 0.0119i 0.491 - 0.526i 



3 LFV processes and muon g — 2 

Even though the soft SUSY breaking parameters are flavor blind at the scale of the SUSY 
breaking mediation, the LFV interactions in the model can induce the LFV source at 
low energies through the renormalization effects. In the following analysis, we assume 
the mSUGRA scenario [19J as the SUSY breaking mediation mechanism. At the original 
scale of the SUSY breaking mediation, we impose the following boundary conditions 
characterized by only five parameters (all of them are assumed to be real in this paper), 
m , M1/2, A , B and \i. Here, m is the universal scalar mass, Mi/ 2 is the universal gaugino 
mass, and A$ is the universal coefficient of the trilinear couplings. The parameters in the 
Higgs potential, B and /i, are determined at the electroweak scale so that the Higgs 
doublets obtain the correct electroweak symmetry breaking VEVs through the radiative 
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breaking scenario [20]. The soft SUSY breaking parameters at low energies are obtained 
through the RGE evolutions with the boundary conditions. 

Although the SUSY breaking mediation scale is normally taken to be the (reduced) 
Planck scale or the string scale (~ 10 18 GeV), in the following calculations we impose 
the boundary conditions at the GUT scale, and analyze the RGE evolutions from the 
GUT scale to the electroweak scale. This ansatz is the same as the one in the so-called 
constrained MSSM (CMSSM). Note that, even with only particle contents introduced in 
the previous section, the minimal SUSY SO(10) model is not the asymptotic free gauge 
theory. When a complete set of the Higgs multiplets is taken into account, the gauge 
coupling would quickly blow up just above the GUT scale. If this is the case, field theory 
description of the model will be no longer meaningful above the GUT scale, and thus it is 
reasonable to introduce the GUT scale cutoff into the model. If there exists new physics 
like M-theory [21], where the string scale comes down to the GUT scale, this treatment 
would be justified. 

The effective theory which we analyze below the GUT scale is the MSSM with the 
right-handed neutrinos. The superpotential in the leptonic sector is given by 

W Y = YpyStiH* + YfieMHi + \m Rij {v R ) t {v R ) 3 + fxH d H u , (8) 

where the indeces i, j run over three generations, H u and H d denote the up- type and down- 
type MSSM Higgs doublets, respectively, and M R . . is the heavy right-handed Majorana 
neutrino mass matrix. We work in the basis where the charged-lepton Yukawa matrix Y e 
and the mass matrix M Ri . are real-positive and diagonal matrices: Y^ = Y ei Sij and M Ri . = 
diag(M Rl , Mr 2 , M Rz ). Thus the LFV is originated from the off-diagonal components of 
the neutrino Dirac Yukawa coupling matrix Y v . The soft SUSY breaking terms in the 
leptonic sector is described as 

" £ soft = 4 ( m !) y h + 4 ( m l) y m + 4 (mf),, e Rj 

+ m 2 Hu HlH u + m 2 H H d H d + (BfiH d H u + \.B V M Ri . v Ri v R j + h.c. 



+ (ASvlJjH* + Ai j e Ri £jH d + h.c. 



+ 1-MxBB + -M 2 W a W a + -M 3 G a G a + h.c) . (9) 
V 2 2 2 y 

As discussed above, we impose the universal boundary conditions at the GUT scale 
such that 

H)y = H)y = K)y = "fai , 

m H u = m H d = m l , 

A^ = AoY^, A^=AoY:\ 

M 1 = M 2 = M 3 = Mi /2 , (10) 
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and extrapolate the soft SUSY breaking parameters to the electroweak scale according 
to their RGEs. The LFV sources in the soft SUSY breaking parameters such as the off- 
diagonal components of (jnfj and A l J are induced by the LFV interactions through the 

neutrino Dirac Yukawa couplings. For example, the LFV effect most directly emerges in 
the left-handed slepton mass matrix through the RGEs such as 



d ( 2 \ d r o 



MSSM 



+ yj^ (mjYX + YXmj + 2Y^mlY v + 2m 2 H YX + 2 4A,) .. (H) 

where the first term in the right hand side denotes the normal MSSM term with no LFV. 
In the leading- logarithmic approximation, the off-diagonal components (z 7^ j) of the 
left-handed slepton mass matrix are estimated as 



where the distinct thresholds of the right-handed Majorana neutrinos are taken into ac- 
count by the matrix L = log[Ma / Mr^S^ . 

The effective Lagrangian relevant for the LFV processes (£i — > £jj) and the muon g — 2 
is described as 

£ e ff = —mutj^F^ (a{P l + A'iP R ) U , (13) 

where Pr,l = (1± 7s)/2 is the chirality projection operator, and Al,r are the photon- 
penguin couplings of 1-loop diagrams in which chargino-sneutrino and neutralino-charged 
slepton are running. The explicit formulas of A^r etc. used in our analysis are summa- 
rized in [22] |2*3*] . The rate of the LFV decay of charged-leptons is given by 

m - Itf) = ^mj, (|Aif + |^| 2 ) , (14) 

while the real diagonal components of Al,r contribute to the anomalous magnetic mo- 
ments of the charged-leptons such as 

5a$> SY = ^ = -ml Re [Af + A%] . (15) 

In order to clarify the parameter dependence of the decay amplitude, we give here an 
approximate formula of the LFV decay rate [22] , 



-.2 



(Amj) 



1 2 



m - try) ~ — m{, x tan 2 /? , (16) 



e m 5 ft 2 ■ v «/y , _ 2 

167T * 167T 2 M| 
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where Ms is the average slepton mass at the electroweak scale, and ( Am|J is the slepton 

mass estimated in Eq. (|12|). We can see that the neutrino Dirac Yukawa coupling matrix 
plays the crucial role in calculations of the LFV processes. We use the neutrino Dirac 
Yukawa coupling matrix of Eq. (JJJ) in our numerical calculations. 

Now we present our results for the rate of the LFV processes and the muon g — 2. 
In Fig. 1, the branching ratio of the muon LFV decay, Br(/i — > ey), is plotted as a 
function of the universal scalar mass m for the various universal gaugino masses, Mi/ 2 = 
400, 600, 800, 1000 GeV, and A = 0. The horizontal line denotes the present experimental 
upper bound, Br(/x — > ej) < 1.2 x 10 -11 [21]. We always take [i > 0, which leads to the 
SUSY contribution to the muon g — 2 being positive. On the other hand, for various 
m = 400, 600, 800, 1000 GeV the branching ratio as a function of Mi/ 2 is depicted in 
Fig. 2. The predicted branching ratios are large so that a parameter space with small 
m and M 1 / 2 has already been excluded. Note that even for m , Mi/ 2 ~ 1 TeV the 
branching ratios well exceed the future planed upper bound, Br(yU — > cy) < 10~ 14 [25 . 
The LFV process may be observed in the near future. The branching ratios for the 
process r — ► /ry are plotted in Fig. 3 and 4 for the same input parameters in Fig. 1 
and 2, respectively, together with the current upper bound, Br(r — > fij) < 0.6 x 10 -6 
In this case, the resultant branching ratio is too small to be accessible to the future 
sensitivity, Br(r — > yrf) < 10~ 9 for input parameters consistent with the current 
bound on the muon LFV decay (see Fig. 1, 2). As can be understood from Eqs. (|T2|) 
and (JIE1), the branching ratio depends on A . For fixed m = 600 GeV and Mi/ 2 = 800 
GeV, the branching ratios of both /i — > ey and r — > /ry are plotted in Fig. 5 as functions 
of Aq. The other LFV processes such as // — > 3e and /i — e conversion rate in nuclei 
are also interesting, since they are strongly constrained by experiments. These processes 
consist of the 1-loop penguin-type and the box-type diagrams. In our parameter space, 
contributions from penguin-type diagrams involving the above Al,r dominates, and we 
can find the approximation formulas, Br(/i — > 3e)/Br(/x — > ej) ~ 7 x 10 -3 [22 j and 
R(fi -> e;Ti(Al))/Br(> -> e 7 ) ~ 5(3) x 10~ 3 (see, for example, Ref. |Zj for detailed 
calculations) for the \i—e conversion with the Ti (Al) target. These are close to the present 
experimental upper bounds, Br(/x — > 3e) < 1.0 x 10 -12 [24] and R(fi —> e; N) < 6.1 x 10 -13 
[2H], when Br(/i — >■ ej) is close to the present upper bound as seen in Fig. 1 and 2. 

Next let us discuss the muon g — 2. Numerical results are depicted in Fig. 6 and 7 for 
the same input parameters in Fig. 1 and 2, respectively. Although these results are almost 
insensitive to the neutrino Dirac Yukawa couplings and the universal trilinear coupling 
Aq, they are correlated with the LFV branching ratios through the dependences on m 
and Mi/2. Note that the input parameters providing Br(/x — > ej) close to the present 
upper bound predict the suitable magnitude for the observed muon g — 2 in Eq. (JTJ). 

If the R-parity is conserved in SUSY models, the lightest superpartner (LSP) is stable. 
The neutralino, if it is the LSP, is the plausible candidate for the cold dark matter (CDM) 
in the present universe. The recent Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite data [2^ provide estimations of various cosmological parameters with greater 
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accuracy. The current density of the universe is composed of about 73% of dark energy 
and 27% of matter. Most of the matter density is in the form of the CDM, and its density 
is estimated to be (in 2a range) 

n CD Mh 2 = 0.11261°;°^ . (17) 

The parameter space of the CMSSM which allows the neutralino relic density suitable for 
the cold dark matter has been recently re-examined in the light of the WMAP data [30 . 
It has been shown that the resultant parameter space is dramatically reduced into the 
narrow stripe due to the great accuracy of the WMAP data. It is interesting to combine 
this result with our analysis of the LFV processes and the muon g — 2. In the case relevant 
for our analysis, tan/3 = 45, /i > and A = 0, we can read off the approximate relation 
between m and My 2 such as (see Figure 1 in the second paper of Ref. |3H] ) 

^o(GeV) = 4^i/ 2 (GeV) + 150(GeV) , (18) 
z8 

along which the neutralino CDM is realized. 4 Mi/ 2 parameter space is constrained within 
the range 300GeV < M 1 / 2 < lOOOGeV due to the experimental bound on the SUSY 
contribution to the b — > 57 branching ratio and the unwanted stau LSP parameter region. 
We show Br(/i — > cy) and the muon g — 2 as functions of M\i 2 in Fig. 8 and 9, respectively, 
along the neutralino CDM condition of Eq. (fTH|) . We find the parameter region, 560GeV < 
^1/2 < 800GeV, being consistent with all the experimental data. 

Before closing this section, we give a comment on the electric dipole moments (EDMs) 
of the charged- leptons. If the diagonal components of A^r have imaginary parts, CP is 
violated and the EDMs of the charged-leptons are given by 

dtje = - mti lm [Af, - A%] . (19) 

These complex A^ r are induced through the renormalization effects in the same manner as 
for the LFV processes. Here the primary source of the CP- violation is the CP-phases in the 
complex neutrino Dirac Yukawa coupling matrix of Eq. (|ZJ) • For the same input parameter 
region analyzed above, we obtain the results of the electron and muon EDMs such as (we 
find d e , < 0) \d e \ = 10~ 34 - 1CT 33 [e cm] and \d^\ = 1(T 31 - 1CT 30 [e cm], respectively, 
which are far below the present experimental upper bounds [21], d e < 1.6 x 10~ 27 [e 
cm] and d M < 3.7 x 10~ 19 [e cm]. As an example of our results, we present the electron 
EDM in Fig. 10 as a function of M\j 2 along the cosmological constraint of Eq. p8jl . If 
the measurement of the electron EDM could be improved by six orders of magnitude, 
d e < 10~ 33 [e cm], as proposed in [HI], there might be a chance to detect the electron 
EDM. 



4 We can see that our resultant soft SUSY breaking mass spectrum at the electroweak scale is almost 
the same as those in the CMSSM. Thus we use the CMSSM result of Eq. (HHJ) as the neutralino CDM 
constraint. 
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4 Summary 



The evidence of the neutrino masses and flavor mixings implies the non-conservation of 
the lepton flavor symmetry in each generations. Thus the LFV processes in the charged- 
lepton sector are allowed. In SUSY theory, there is the possibility that the rate of the LFV 
processes will be accessible to the current or future experimental bounds. The concrete 
information of the LFV interactions is necessary to evaluate the rate of the LFV processes. 

It has recently been shown [3] that the minimal SUSY SO(10) model can simultane- 
ously accommodate all the observed quark-lepton mass matrix data involving the neutrino 
oscillation data with appropriately fixed free parameters. In this model, the neutrino 
Yukawa coupling matrix are completely determined, whose off-diagonal components are 
the primary source of the lepton flavor violation in the basis where the charged-lepton 
and the right-handed neutrino mass matrices are real and diagonal. Using this Yukawa 
coupling matrix, we have calculated the rate of the LFV processes assuming the mSUGRA 
scenario. The branching ratio of the \x — > ej process is found to be large so that a pa- 
rameter space with small universal scalar and gaugino masses has been already excluded 
by the present experimental upper bound. Even for 1 TeV input parameters, we found 
that the branching ratio well exceeds the future planed experimental upper bound. 

We also have calculated the SUSY contribution to the muon g — 2, which is correlated 
with the rate of the LFV processes through the input parameters in mSUGRA scenario. 
The resultant magnitude of the muon g — 2 is found to be within the range of the recent 
result of Brookhaven E821 experiment for the input parameters providing the \x — > ey 
branching ratio close to the current experimental bound. 

In CMSSM the parameter region realizing the neutralino dark matter scenario has 
been dramatically reduced by the recent WMAP data. We have performed the same 
analysis for the LFV processes and the muon g — 2 by restricting the input parameters 
within the cosmologically allowed region. We have found that there exists a parameter 
region consistent with all the data. 
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Figure 1: The branching ratio, Log 10 [Br(/i — > e^y)], as a function of mo (GeV) for Mi/ 2 = 
400, 600, 800, 1000 GeV (from top to bottom) with A = and fj, > 0. 




Figure 2: The branching ratio, Log 10 [Br(/i — > ej)], as a function of Mi/ 2 (GeV) for 
m = 400, 600, 800, 1000 GeV (from top to bottom) with A = and /i > 0. 
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Figure 3: The branching ratio, Log 10 [Br(r — > /xy)], as a function of mo (GeV) for Mi/ 2 = 
400, 600, 800, 1000 GeV (from top to bottom) with A = and y. > 0. 



-5 




400 600 800 1000 1200 
Mi/2 [GeV] 

Figure 4: The branching ratio, Log 10 [Br(r — >■ 7x7)], as a function of Mi/ 2 (GeV) for 
m = 400, 600, 800, 1000 GeV (from top to bottom) with A = and y > 0. 
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Figure 5: The branching ratios, Log 10 [Br(r — > fxj)] (top) and Log 10 [Br(/i — > ej)] (bot- 
tom) as functions of A (GeV) for m = 600 GeV and M 1/2 = 800 GeV. 
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Figure 6: The SUSY contribution to the muon g — 2 in units of 10 as a function of mo 
(GeV) for M 1/2 = 400, 600, 800, 1000 GeV (from top to bottom) with A = and fi > 0. 
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Figure 7: The SUSY contribution to the muon g — 2 in units of 10~ 10 as a function of 
Mi/2 (GeV) for m = 400, 600, 800, 1000 GeV (from top to bottom) with A = and 
H > 0. 
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Figure 8: The branching ratio, Log 10 [Br(/x — > e7)], as a function of M1/2 (GeV) along 
the cosmological constraint of Eq. ()18|) . 
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Figure 9: The SUSY contribution to the muon g — 2 in units of 10~ 10 as a function of 
Mi/2 (GeV) along the cosmo logical constraint of Eq. ([180. 
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Figure 10: The electron EDM, Log 10 [|rf e |[e cm]], as a function of Mi/ 2 (GeV) along the 
cosmo logical constraint of Eq. (fTS[) . 
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